Localization-delocalization transition in 2D quantum percolation model 



Md Fhokrul Islam and Hisao Nakanishi 
Department of Physics, Purdue University, West Lafayette, IN47907 
(Dated: February 2, 2008) 

We study the hopping transport of a quantum particle through randomly diluted percolation 
clusters in two dimensions realized both on the square and triangular lattices. We investigate the 
nature of localization of the particle by calculating the transmission coefficient as a function of 
energy (—2 < E < 2 in units of the hopping integral in the tight-binding Hamiltonian) and disorder, 
q (probability that a given site of the lattice is not available to the particle). Our study based on 
finite size scaling suggests the existence of delocalized states that depends on energy and amount 
of disorder present in the system. For energies away from the band center (E = 0), delocalized 
states appear only at low disorder (q < 15%). The transmission near the band center is generally 
very small for any amount of disorder and therefore makes it difficult to locate the transition to 
delocalized states if any, but our study does indicate a behavior that is weaker than power-law 
localization. Apart from this localization-delocalization transition, we also find the existence of two 
different kinds of localization regimes depending on energy and amount of disorder. For a given 
energy, states are exponentially localized for sufficiently high disorder. As the disorder decreases, 
states first show power-law localization before showing a delocalized behavior. 



I. INTRODUCTION 

Unlike a classical particle, the transport of a quan- 
tum particle through a system is greatly influenced by 
the interference and tunneling effects. While the avail- 
ability of a spanning path is the sole criteria for the 
transmittance of a classical particle through a system, a 
quantum particle may exhibit zero or very low transmit- 
tance even for a completely ordered system depending 
on details such as the boundary condition or the energy 
of the particle^. The quantum percolation system that 
we have investigated here includes the interference ef- 
fect but does not include the tunneling effect. We, thus, 
expect a higher connectivity in underlying geometry to 
be required for non-zero transmission compared to its 
classical counterpart. 

A major motivation for studying such a system is 
the question of whether a localizcd-to-dclocalized (or 
perhaps, metal-to-insulator) transition exists in a two- 
dimensional (2D) system. The Anderson model and 
the quantum percolation model are two of the more 
common theoretical models that are used to study the 
transport properties of disordered systems. While the 
literature on both models agree on the existence of such 
a transition in three dimensions'^, the same question 
for quantum percolation in two dimensions appears to 
have remained a subject of controversy for over two 
decades. Based on the one-parameter scaling theory 
of Abrahams et al£, it was widely believed that there 
can be no metal-to-insulator transition in 2D univer- 
sally in the absence of a magnetic field or interactions 
for any amount of disorder. Moreover, the scaling the- 
ory predicts that all states are exponentially localized 
in thermodynamic limit for any amount of disorder and 
therefore that no transition exists except at zero dilu- 
tion. (However, see Goldenfeld and Haydock 6 which 
asserts the existence of a transition between two dif- 
ferent kinds of localized regimes at a finite disorder in 
addition to a localizcd-delocalized transition at an in- 



finitesimal dilution, even for the Anderson model in two 
dimensions.) 

However, whether the scaling theory also applies 
to quantum percolation has been debated in recent 
years. Even restricting attention to quantum percola- 
tion which lacks many effects that are expected to play 
important roles in metal-insulator transitions, there is 
a long-standing controversy as to the presence or ab- 
sence of an extended state and of a phase transition be- 
tween the prevalent localized state and a more elusive 
extended state in two-dimensions. On one hand, some 
studies such as those made using the dlog Pade approx- 
imation method^, real space renormalization method 8 , 
and the inverse participation ratio^ found a transition 
from exponentially localized states to non-exponentially 
localized states for a range of site concentrations be- 
tween 0.73 < p q < 0.87 on the square lattice. So did 
a study of energy level statistics^, one of the spread 
of a wave packet initially localized at a siteii, and one 
of a transfer matrix 1 -, where the nature of the delocal- 
ized state remained not fully understood. On the other 
hand, studies such as the scaling work based on numer- 
ical calculation of the conductance^, the investigation 
of vibration-diffusion analogy 14 , finite-size scaling anal- 
ysis and transfer matrix methods 1 ^, and vector recur- 
sion technique 16 found no evidence of a transition. A 
study by Inui et at 17 found all states to be localized 
except for those with particle energies at the middle of 
the band and when the underlying lattice is bipartite, 
such as a square lattice. More recently, Cuansing and 
Nakanishi 7 used an approach first suggested by Daboul 
at al.' to calculate conductance directly for clusters of 
up to several hundred sites and, extrapolating those re- 
sults by finite-size scaling, suggested that delocalized 
states exist and thus a transition would have to exist as 
well. The current paper extends the relevant portion of 
the latter work by studying much larger clusters, which 
has allowed considerably more detailed analyses. 

In the mean time, experiments performed in early 



1980's on different 2D system a 19 ' 20 ' 21 confirmed the 
scaling theory predictions. However, a number of ex- 
periments on dilute low disordered Si MOSFET and 
GaAs/AlGaAs heterostructures that appeared more re- 
cently seem to suggest that a metallic state may be 
possible in two dimension s 22 ! 23 ' 24 ' 2 ^. For reviews of 
these experiments, see Abrahams et al£L and references 
therein. 

In this work, we do not address the issues of these ex- 
periments, but rather concentrate on the formally much 
simpler quantum percolation model which has neither 
magnetic field nor interactions but contains binary dis- 
order with infinite barriers at randomly diluted sites. 
Previously we have investigated the same problem using 
a dynamical approach where we have studied the prop- 
erties of a disordered system by tracking how a quan- 
tum particle, described by a wave packet, propagates 
through the system 2 ^. In this paper we adopt a station- 
ary state approach where we calculate the transmission 
characteristics by solving time independent Schrodinger 
equation for tight binding Hamiltonian. In Section II we 
describe the model and the numerical approach used in 
this work. In Section III we discuss the numerical re- 
sults and in Section IV we present the summary and 
conclusion of our study. 



only one lattice site on the opposite edge of the clus? 
ter. Another possible connection type is the busbar 
type contact, where all the lattice points on the input 
side edge of the cluster are connected to the input lead, 
while all the lattice points on the output side of the 
cluster are connected to the output lead. Figs. [1] and [2] 
illustrate the connection of the leads for the square and 
the triangular lattices respectively. 
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II. QUANTUM PERCOLATION MODEL AND 
NUMERICAL APPROACH 

We study quantum percolation that is described by 
the Hamiltonian 

H= J2 Vij\i)U\+h.c (1) 

<ij> 

where i) and \j) are tight binding basis functions 
at sites i and j, respectively, and Vij is the hopping 
matrix element which is equal to zero if i and j are 
not nearest neighbors. We have realized this model on 
both square and triangular lattices that can have at 
most 4 and 6 nearest neighbors, respectively. If the 
system is completely ordered, Vij = Vq (uniform) and 
Vb sets the overall energy scale, where we use Vq = 1 as 
the nominal standard value. On the other hand, since 
in this work we are interested in transport through a 
disordered system, we will introduce random dilution 
by removing a fraction of sites from the lattice and set 
Vij = for the bonds between the diluted sites and their 
neighbors. Vij = for nearest neighbors i and j means 
that an infinite barrier exists between the pair of sites. 

To study the transmission of a quantum particle we 
connect two semi-infinite ID leads, one as the input 
and the other as the output lead, to the 2D cluster. Al- 
though the system can be studied using different types 
of connection of the leads, in this study we only use a 
point-to-point type contact where the input lead is con- 
nected to only one lattice site on the input side edge 
of the cluster and the output lead is connected also to 



(a) (b) 

FIG. 1: 3 x 3 Square Lattice: (a) point-to-point connection 
and (b) busbar type connection. The letters label the lattice 
points of the cluster part of the Hamiltonian, while numbers 
label those of the leads. The same sequence of labeling is 
used for all sizes of the clusters in this work. 
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(a) (b) 

FIG. 2: 3 x 3 Triangular Lattice: (a) point-to-point connec- 
tion and (b) busbar type connection 

The wave function of the entire cluster-lead sys- 
tem can be calculated by solving the time independent 
Schrodinger equation: 



3 



Hip = Ei/; 

(Win \ 
^cluster 
Ipout ) 

and Vm = {^-(n+i)} and ip out = {ip+(n+i)}, 
n = 0, 1, 2 . . ., are the input and output chain part of 
the wave function respectively. 

Since the leads are of infinite length, the matrix form 
of the Schrodinger equation (Eq. [2|) becomes an infinite 
size problem. To reduce it to a numerically finite prob- 
lem we use an ansatz proposed by Daboul et alX which 
assumes that the input and output part of the wave 
function are of the form of plane waves: 



VW^+(n+l) = fe"" [ ' 

where r is the amplitude of reflected wave, t is the 
amplitude of the transmitted wave. This ansatz is con- 
sistent with the Schrodinger equation only for the wave 
vector q that is related to the energy E by 

E = e~ iq + e lq (4) 

Using this ansatz along with the energy restriction 
Eq. Q , the matrix equation for a 3 x 3 cluster connected 
to semi-infinite chains (Fig. QJi) reduces to (for details 
see reference^) 
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Here c is the coupling of the leads with the cluster 
and its value is set to 1 for all the calculations done 
in this work. The effect of c ^ 1 on transmission is 
discussed in one of our previous works2£. The busbar 
configuration of Fig. [1] (b) has a similar expression. 

The Eq. ([5]) is the exact expression for a 2D system 
connected to semi-infinite chains with continuous eigen- 
values ranging between -2 and +2. The spectrum is 
continuous because it is still effectively infinite and it is 
non-degenerate except for the reversal of left and right. 

The main advantage of using this ansatz is that it 
not only allows us to calculate the wave function but 
also helps us to study the transmission characteristics 
of the corresponding state directly. The transmission 
and the reflection coefficients are obtained by taking 
the absolute square of t and r respectively, ie T = \t\ 2 
and R = \r\ 2 . A disadvantage on the other hand is 
that Eq. (|1]) that relates the wave vector of the incident 
particle with the energy, restricts the energy of the par- 
ticle to between -2 and +2. This restricts our ability 
to study the system in the whole possible energy range 
since for our system the energy could in principle range 
between -4 and +4. This is due, of course, to the effec- 
tively one-dimensional nature of the system forced by 



attaching ID semi-infinite leads and by looking at plane 
waves spreading over the entire leads. 



III. NUMERICAL RESULTS 

The model that we are using to calculate transmis- 
sion has two adjustable parameters, namely the energy 
of the particle, E and the amount of disorder present in 
the system, q (probability that a given site is not avail- 
able to the hopping particle). To study the presence or 
absence of a localization-delocalization transition, one 
needs to investigate the behavior of the 2D system in 
the thermodynamic limit. This, however, is not possi- 
ble in numerical methods, and therefore, we resort to 
a finite size scaling approach in which we calculate the 
transmission while gradually increasing the size of the 
system for a given amount of disorder. The result is 
then extrapolated to study the bulk behavior in the 
thermodynamic limit. 

Though there are many combinations possible for 
connecting the input and the output leads with the 
cluster, for the work that follows, the leads are con- 
nected diagonally with the clusters, since this arrange- 



ment generally gives higher transmission. In addition 
to that, to minimize the effect of the boundary on the 
interior property of the disordered clusters, we made 
good contacts by keeping the nine sites nearest to both 
the input and the output contact points always occu- 
pied (that is, available to the hopping particle). We 
calculate the transmission as a function of the system 
size for many levels of disorder and at different energies. 
The general trend for all the transmission curves at dif- 
ferent energies are similar, so we will discuss here only 
two energies, one that is away from the band center and 
one very close to the band center. 



have calculated the transmission by gradually increasing 
the size of the clusters from 10 x 10 to a maximum size of 
180 x 180. We have randomly generated one thousand 
clusters of a given size for each level of disorder and 
average transmission is calculated for each size of the 
cluster and each level of disorder. The log-log plot of 
transmission against the size of the clusters at E = 1.6 
is shown for various disorder levels in the Fig. [3] 



A. Energies away from the band center 

We first study the transmission at energy E = 1.6 for 
different levels of disorder. For each level of disorder, we 




FIG. 3: The log-log plot of transmission through disordered clusters at E — 1.6. Each data point is the average over 1000 
independent realizations. The transmission curves for 5%, 10%, 13%, 15% and 20% disorders are separately shown in the 
inset. 



We observe from Fig. [3] that at lower disorder trans- 
mission decreases almost linearly in the log-log plot as 
the size of the clusters increases, suggesting a power-law 
behavior of transmission in linear scale. This trend con- 
tinues until disorder, q, increases to about 28%. Above 
that, transmission falls exponentially as is evident from 
the lowest three curves in Fig. [3j We have fitted the 
data both to power-laws and exponentials and the best 
fit forms along with the corresponding fitting param- 
eters and coefficient of regression, i?| 2 , for different 
levels of disorder are presented in Table |TJ Although 
there are many data points that appear to be several 
cr's from the best fit exponentials particularly for higher 
disorder, they are believed to be due to the discrete 
and loose-packed structure of the lattice and not due 
to some unknown systematic errors. (The transmission 
curves of the 13% and 26% disorder are not shown in 
the Fig. [3] since they are very close to the nearby trans- 
mission curves). 



TABLE I: Table for fitting parameters at E = 1.6. Shown 
in the parentheses are the lower and upper bounds for 95% 
confidence level. 



q 


Fit 
equation 


Parameters 


|i?| 2 


a 


b 


5% 


T = a-L~ b 


0.78 (0.69, 0.89) 


0.26 (0.23, 0.30) 


0.95 


10% 


1.09 (0.89, 1.32) 


0.44 (0.39, 0.50) 


0.96 


13% 


1.24 (1.01, 1.52) 


0.54 (0.49, 0.60) 


0.97 


15% 


1.404 (1.16, 1.70) 


0.63 (0.58, 0.68) 


0.98 


20% 


2.70 (2.25, 3.25) 


0.96 (0.91, 1.01) 


0.99 


23% 


4.08 (3.47, 4.81) 


1.19 (1.15, 1.23) 


0.99 


25% 


11.3 (6.20, 20.8) 


1.59 (1.43, 1.75) 


0.97 


26% 


16.3 (8.98, 29.5) 


1.77 (1.61, 1.93) 


0.97 


28% 


36.2 (18.4, 71.5) 


2.12 (1.93, 2.30) 


0.98 


30% 


T = a- e~ bL 


0.18 (0.12, 0.25) 


0.06 (0.05, 0.06) 


0.97 


35% 


0.24 (0.12, 0.46) 


0.12 (0.11, 0.13) 


0.97 


38% 


0.96 (0.15, 6.12) 


0.22 (0.19, 0.26) 


0.94 



It should be noted that the curve fitting is performed 
in log-log scale and then converted to linear scale along 
with the parameters. Thus to fit the data, say for 
q = 15%, we first calculate the logarithm of the size 
of the cluster and the corresponding transmission, ic 
x = log(L) and y = log(T). We then fit these logarith- 
mic data with the curve, y — m.x + c, and converted 
the result back to linear scale. 

Because of the large differences in transmission for 
different disorder, we used the logarithmic scale for 
transmission in Fig. [3j which in turn makes it difficult 
to see the finer details of the transmission curve partic- 
ularly at lower disorders. For instance, Fig. [3J appears 
to suggest an excellent linear fit at disorders below 15%. 
To investigate more closely the nature of transmission 
at lower disorder we have plotted the first five curves 
of Table Q] separately in the inset of Fig. [3J It is evi- 
dent from the inset that transmission decreases much 
more slowly than a power-law at larger values of L for 



disorder up to about 13%. Above 13% power-laws do 
appear to give good fits of the data for the entire range 
of cluster sizes used. To investigate this lower disorder 
regime further, we have obtained the transmission for 
2% and 3% disorder and have plotted these data along 
with those for 5% to 15% disorder in a linear graph as 
shown in Fig. Q] for larger clusters (L > 25). Each data 
set is fitted with two different curves namely, 
a power-law: 

T = a-L- b (6) 
and an exponential with offset: 

T = a-e- bL + c (7) 
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FIG. 4: The linear plot of transmission through clusters with 
2%, 3%, 5%, 10%, 13% and 15% disorders at E = 1.6 for 
larger sizes of the clusters (L > 25). Each data point is 
the average over 1000 independent realizations. The dotted 
lines represent power-law fits and the solid lines represent 
the exponential fits with offset. 

The dotted and solid lines represent fits with Eq. ([6]) 
and (J7J), respectively. As we can see from Fig. HI ex- 
ponential with offset, c, gives significantly better visual 
fits for the data compared to power-laws at least for 2% 
and 3% disorder, while the goodness of fits for higher 
disorders appear to be a toss-up. 

To compare the goodness of fits we need to consider 
not only the visual fits but the values of |i?| 2 and SSE 
(sum square error) as well, and it is evident from Ta- 
ble [TT] and Fig. 2b that an exponential with offset gives 
significantly better fits than a power law, except for 15% 
disorder, at which both fits are quite close. 

The above study, thus, suggests the existence of three 
regimes of transmission. At higher disorder transmis- 
sion drops exponentially as the size of clusters increases. 



TABLE II: Table for fitting parameters of the linear plot at 
E — 1.6. Shown in the parentheses are the lower and upper 
bounds for 95% confidence level. 



q 


T=a-L~ b 


T = a- e~ bL + c 


a 


b 


1 D\2 

\R\ 


SSE 


a 


b 


c 


I r>l2 

\R\ 


SSE 


noy 
2% 


0.78 


0.22 


0.79 


5x 
10~ 3 


0.48 


0.05 


0.27 
(0.26, 0.29) 


0.89 


2x 
10 -3 


ooy 

3 To 


0.67 


0.20 


0.91 


lx 

i n-3 
10 


0.27 


0.03 


0.25 
(0.24, 0.26) 


0.96 


5x 

10 


07o 


0.61 


0.21 


0.96 


3x 
10" 4 


0.17 


0.02 


0.21 
(0.19, 0.23) 


0.95 


4x 
10~ 4 


10% 


0.78 


0.37 


0.97 


3x 
10 -4 


0.21 


0.02 


0.12 
(0.11, 0.14) 


0.97 


3x 
10" 4 


13% 


0.95 


0.48 


0.96 


5x 
10" 4 


0.22 


0.03 


0.09 
(0.06, 0.1) 


0.94 


6x 
10" 4 


15% 


1.00 


0.55 


0.96 


4x 
10" 4 


0.18 


0.02 


0.05 
(0.02, 0.08) 


0.96 


4x 
10" 4 



At intermediate disorder, transmission follows powe*? 
law, whereas at low disorder transmission approaches 
a constant offset, c, suggesting possible delocalization 
of the states. The value of the offset, however, depends 
on energy. For a given disorder, c decreases as the en- 
ergy increases, except for energies very close to the band 
center where transmission itself is very low. 



B. Energies near the band center 

In this subsection we discuss the nature of transmis- 
sion at an energy near the band center. We investigate 
the scaling behavior at E — 0.001 for different levels of 
disorder. The log-log plot of transmission as a function 
of system size is shown in the Fig. [5J 




FIG. 5: The log-log plot of transmission through disordered clusters at E = 0.001. Each data point is the average over 1000 
independent realizations. The transmission curves for 5%, 6%, 7%, and 8% disorders are also separately shown in the inset. 



At this energy, transmission decreases exponentially 
for disorder as low as 9% as is evident from Fig. [5l 
Below 9% the transmission curves appear to fit well 
with power-laws. The fitting parameters of the different 
curves are tabulated in the Table Ifffl along with the 
coefficient of regression, |i?| 2 . 



TABLE III: Table for fitting parameters at E = 0.001. 
Shown in the parentheses are the lower and upper bounds 
for 95% confidence level. 



q 


Fit 
equation 


Parameters 


\R\ 2 


a 


b 


5% 


T = a-L- b 


38.2 (22.7, 64.5) 


1.43 (1.29, 1.57) 


0.97 


6% 


70.0 (69.9, 113) 


1.70 (1.57, 1.83) 


0.98 


7% 


204 (97.6, 429) 


2.12 (1.93, 2.32) 


0.97 


8% 


691 (216, 2206) 


2.62 (2.31, 2.93) 


0.96 


9% 


T = a- e~ bL 


0.85 (0.53, 1.35) 


0.06 (0.06, 0.07) 


0.96 


10% 


1.88 (1.31, 2.69) 


0.10 (0.11, 0.10) 


0.99 


15% 


8.00 (1.90, 33.7) 


0.24 (0.21, 0.26) 


0.97 



The |i?| 2 values for 5%, 6%, 7% and 8% disorder from 
the Table Irnl suggest a good fit of these data with power- 
laws. However, if we take a closer look at the transmis- 
sion curves for these disorders in the inset, we observe 
a significant deviation from straight lines as the size of 
the cluster increases. It is clear that at lower disorders, 
the falling trend in transmission is much slower than the 
power-law similarly to what we have observed at ener- 
gies far from the band center. However, transmission is 
generally much smaller for the same amount of disorder 
compared with the case of energies away from the band 
center. 

As before we have obtained two more transmission 
data for 2% and 3% disorder and the results are plotted 
in Fig. [SJi along with that for 5% disorder (the trans- 
mission curve for 6% disorder is too close to that of 5% 
disorder and is not shown in the the figure). We show in 
Fig. [5^ transmission data for the entire range of cluster 
sizes. We have fitted these data with both a power-law 
and an exponential with offset as before, and the best 
fit curves along with the fit parameters are tabulated 
in the Table ITVl The dotted lines represent the power- 
law fits, and the solid lines represent exponential fits 
with offset. We notice that at this energy the best fit 
curve for each of the disorder amounts shown is clearly 
an exponential with offset rather than a power-law. For 
more direct comparison with Fig. U where we consid- 
ered only the larger clusters (L > 25), we also show 
the corresponding figure and fits for the same range of 
cluster sizes in Fig. m> as well. In the latter figure, the 
difference in goodness of fits between power-law fits and 
exponential fits with offset is not as large, but it is still 
clear that the latter fits the data better. 

The issue of localization length, L;, must also be dis- 
cussed. Our finite size scaling approach does not rely on 
any single system size but rather focuses on the trend as 
it increases. In fact this analysis detects different trends 
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FIG. 6: The linear plot of transmission through clusters 
with 2%, 3% and 5% disorder at E = 0.001. (a) plot with 
all data points and (b) plot with data points only for larger 
sizes of the clusters. Each data point is the average over 
1000 independent realizations. The dotted lines represent 
the power-law fits. 

depending on the dilution amount and we can estimate 
the localization length from the data themselves at least 
in the clearly observed exponentially localized regime. 
Table [V] summarizes the simulation estimates of L\ at 
different energies in the regions of disorder where the re- 
spective exponential fits (a-e~ bL ) are significantly better 
than other types of fits. The localization length in the 
exponential regime can be estimated from Li ~ b^ 1 . It 



TABLE IV: Table for fitting parameters of the linear plot 
at E — 0.001. Shown in the parentheses are the lower and 
upper bounds for 95% confidence level. 



q 


T = a- L~ b 


T = a- e~ bL + c 


a 


b 


1 n 1 2 

\R\ 


SSE 


a 


b 


c 


\R\ 


SSE 


2% 


9.25 


0.69 


0.96 


2x 
10 -2 


1.19 


0.02 


0.15 
(0.09, 0.20) 


0.99 


3x 
10 -3 


3% 


22.7 


1.03 


0.99 


7x 
10 -3 


1.44 


0.03 


0.11 
(0.09, 0.12) 


0.99 


lx 
10 -3 


5% 


141 


1.73 


0.99 


2x 
10 -3 


1.79 


0.05 


0.03 
(0.02, 0.04) 


0.99 


5x 
10" 4 


6% 


198 


1.95 


0.99 


4x 
10" 4 


1.53 


0.06 


0.02 
(0.01, 0.03) 


0.99 


7x 
10" 4 



is evident that most of our system sizes are sufficiently 
large compared with these estimates of Li and thus our 
results are internally consistent with exponential local- 
ization at these amounts of disorder. 



TABLE V: Estimated values of the localization lengths from 
T — a ■ e~ bL , localization length being b^ 1 



q 


Localization length (in units of lattice constant 




E = 


E = 


E = 


E = 


E 


E = 


E = 




0.001 


0.05 


0.5 


1.0 


\:> 


1.6 


1.9 


9% 


15.63 














10% 


9.71 














15% 


4.26 














20% 


2.69 


16.67 












25% 




6.29 










27.03 


30% 








19.30 




17.24 


13.26 


31% 






13.89 




16.69 






32% 






13.89 










35% 






8.85 


6.75 


7.12 


8.16 


5.37 


38% 






4.85 


4.32 


4.14 


4.44 





For completeness we present here the scaling result 
for triangular lattice only for energy E — 1 for different 
disorder. The curves appear very similar to those for 
square lattice as is evident from the Fig. [7| except that 
the exponentially localized regime appears at a higher 
disorder. The result is expected since triangular lattice 
has six nearest neighbors and therefore requires more 
disorder to reduce the transmission. 



IV. SUMMARY AND CONCLUSION 

In this paper we have studied the behavior of a quan- 
tum particle in two dimensional disordered clusters in 
quantum percolation model. Our approach is based on 
calculating the transmission of the particle that enters 
into the cluster through a one dimensional lead at one 
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FIG. 7: The log-log plot of transmission through disordered 
clusters in triangular lattice at E = 1. Each data point is 
the average over 1000 independent realizations. 

side of the cluster and exit through another lead at the 
other side of the cluster. 

Our study based on finite size scaling suggests the 
existence of three different regimes depending on disor- 
der. The range of these regimes, however, depends on 
the energy of the particle. Although we showed in pre- 
vious sections results from two representative energies 
(E = 1.6 that is far from the band center and E — 0.001 
that is very close to it), we have in fact obtained and 
analyzed the data for many more values of the energy. 
Table IVT1 shows our approximate estimates of the range 
of disorder for which these regimes exist as the energy 
is varied. 

TABLE VI: Classification of different regime in 2D disor- 
dered system. 



Energy 


Range of disorder 




delocalized 


power law 


exponential 




states 


localization 


localization 




regime 


regime 


regime 


0.001 


0-6% 


7% - 8% 


>9% 


0.05 


0- < 15% 


> 15%- < 20% 


> 20% 


0.5 


- 15% 


> 15% - 31% 


> 32% 


1.0 


- 15% 


> 15% - 28% 


> 30% 


1.2 


- 15% 


> 15% - 30% 


> 31% 


1.6 


0- < 15% 


> 15% - 28% 


> 30% 


1.9 


0- < 15% 


> 15% - 23% 


> 25% 



Our study, thus, suggests that at lower disorder states 
are delocalized, contrary to the one parameter scaling 
theory of Abrahams et a£. In addition to delocalized 



transition we observe a second kind of phase transi- 
tion between the power-law to exponentially localized 
regimes. 

For all energies considered in this work, except for 
those close to E = 0, the delocalized states appear if 
the disorder is less than about 15%. Above 15%, there 
exist two different localization regimes. For interme- 
diate energies (0.5 < E < 1.6), states of the particle 
show a weaker form of localization characterized by a 
power-law size dependence of transmission for the range 
of disorder between 15% to ~ 30%. Above this range all 
states are localized exponentially. For energies close to 
E — 2, the width of this power-law localization regime is 
reduced and states become exponentially localized even 
for 25% disorder. 

The transmission near the band center, however, dif- 
fers in some important ways compared to that at other 
energies. At E ~ 0, transmission is very low in the ther- 



modynamic limit even for very low disorder. Though 1 
it is difficult to locate the derealization transition be- 
cause of the low transmission, our finite size analysis 
does show a behavior much weaker than power-law lo- 
calization (inset of Fig. [5]) . The width of the power-law 
localization regime is very small and exponentially lo- 
calized states appear for disorder as low as 9%. 
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